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Long XMM observation of the Narrow-Line Seyfert 1 galaxy 
IRAS 13224-3809: rapid variability, high spin and a soft lag 
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ABSTRACT 

Results are presented from a 500 ks long XMM-Newton observation of the Narrow-Line 
Seyfert 1 galaxy IRAS 13224-3809. The source is rapidly variable on timescales down to a few 
100 s. The spectrum shows strong broad Fe-K and L emission features which are interpreted as 
arising from reflection from the inner parts of an accretion disc around a rapidly spinning black 
hole. Assuming a power-law emissivity for the reflected flux and that the innermost radius 
corresponds to the innermost stable circular orbit, the black hole spin is measured to be 0.988 
with a statistical precision better than one per cent. Systematic uncertainties are discussed. A 
soft X-ray lag of 100 s confirms this scenario. The bulk of the power-law continuum source is 
located at a radius of 2-3 gravitational radii. 

Key words: accretion, accretion discs - black hole physics - line:profiles - X-rays:general 



1 INTRODUCTION 

The X-ray emission from most Narrow-Line Seyfert 1 galaxies 
(NLS1) is characterised by a steep soft X-ray spectrum and rapid 
variability. The most extreme such objects are 1H 0707-495 and 
IRAS 13228-3809, which both show a sharp drop above 7 keV in 
XMM Newton spectra (Boiler et al 2002, 2003). 1H 0707-495 has 
been further studied several times with XMM, including a long 
500 ks dataset in 2008 which revealed broad iron K and L lines and 
a soft lag of about 30 s (Fabian et al 2009; Zoghbi et al 2010, 201 1). 
In contrast, IRAS 13224-3809 only had 64 ks of XMM data (Ponti 
et al 2010; Gallo et al 2004), despite showing spectacular variabil- 
ity during ROSAT (Boiler et al 1997) and ASCA (Dewangan et al 
2002) observations. New observations totalling 500 ks have now 
been made with XMM in 201 1 and reported here. 

The unusual spectrum and 7 keV drop of both objects have 
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been interpreted as due to either intervening absorption or strong 
relativistic blurring of a reflection component (Boiler et al 2002, 
2003; Fabian et al 2004). 1H 0707-495 dropped into a low state for 
about 2 months at the start of 201 1 during which an XMM spectrum 
showed evidence for even more blurring. The results are consistent 
with the power-law component of the X-ray source lying within one 
gravitational radius of the central black hole (Fabian et al 2012). In 
the normal state, one third of this component extends to ~ 20r g . 

The combination of the above results with the reverbera- 
tion lags in 1H 0707-495 and in over a dozen other sources 
(Emanoulopolous, McHardy & Papdakis 2010; Tripathi et al 201 1; 
De Marco et al 2011, 2012; Zoghbi & Fabian 2011 and Zoghbi et 
al 2012) provides very strong support for the reflection model for 
the X-ray emission of Seyfert galaxies. In this model the primary 
power-law component lies above the inner accretion disc around 
the black hole and produces the X-ray reflection component by ir- 
radiation of the disc (see e.g. Fabian & Ross 2010). The soft lags 
are then the light travel time difference between the power-law and 
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reflection components as detected by the observer. The new data 
presented in this paper are interpreted within the reflection model. 

IRAS 13224-3809 is a radio quiet (1.4 GHz flux of 5.4 mjy, 
Feain et al 2009) NLS1 at redshift z = 0.066. For a flat ACDM 
cosmology with Hq = 71 km s _1 Mpc -1 , its luminosity distance is 
293 Mpc. 



2 OBSERVATIONS AND DATA REDUCTION 

IRAS 13224-3809 was observed for ~ 500 ks with the XMM- 
Newton satellite (Jansen et al 2001) from 2011 July 19 to 
2011 July 29 (Obs. IDs 0673580101, 0673580201, 0673580301, 
0673580401). We focus on the data from the EPIC-pn camera 
(Striider et al 2001). The first observation was taken in full win- 
dow imaging mode, and the following three in large window imag- 
ing mode. All of the data were reduced in the same way, using 
the XMM-Newton Science Analysis System (SAS v. 1 1.0.0) and the 
newest calibration files. 

The data were cleaned for high background flares, resulting in 
a final total exposure of 300 ks. The data were selected using the 
condition PATTERN < 4. Pile-up effects were not significant in any 
of the observations. 

The source spectra were extracted from circular regions of ra- 
dius 35 arcsec, which were centered on the maximum source emis- 
sion, and the background spectra were chosen from a circular re- 
gion of the same size and on the same chip. The position of the 
background regions were chosen to avoid the Cu-K emission lines 
from the electronic circuits behind the pn CCD that contaminate 
the background at 8.0 and 8.9 keV. The response matrices were 
produced using RMFGEN and ARFGEN in SAS. 

The spectra from the four observations were merged before 
fitting using MATHPHA in FTOOLS, and the resulting combined 
spectrum was rebinned to contain a minimum of 20 counts per 
bin. Spectral fitting was performed using XSPEC v 12. 5.0 (Arnaud 
1996). Quoted errors correspond to 90% confidence level. Energies 
are given in the rest frame of the source. Quoted abundances refer 
to the solar abundances in Anders & Grevesse (1989). 



3 LIGHTCURVE AND VARIABILITY 

The lightcurve of the long XMM observation in 2011 is shown in 
Fig. 1. Start and stop dates are 201 1 July 19 and 29. The source is 
clearly highly variable with several pronounced upward spikes of 
emission. The observation consists of 4 orbits of XMM, the first 
two of which are contiguous. Gaskell (2003) used ASCA data to 
show that the X-ray light curves can be lognormal, in the sense that 
a frequency histogram of log(count rate) is gaussian. We find that 
this is only a fair description of Orbit 1 and is a poor description of 
the Orbits 3 and 4. 

The light curve is reminiscent of that seen with the ROSAT 
HRI, reported by Boiler et al (1997). Observations were made ev- 
ery day for a month resulting in a total exposure of ~ 1 10 ks and 
mean observation length ~ 3 ks. 5 flares of emission were found 
rising above about 150 ct ks -1 . Using WEBPIMMS (Mukai 1993) 
and a 0.1(0.2) keV blackbody model (see next Section for detailed 
spectrum) we find that 150 ct ks~' corresponds to about 6(3.4) 
EPIC pn ct s _1 . We see 3 flares above 6 pn ct a- -1 and ~ 11 above 
3.4 pn ct s in 500 ks. The numbers depend on the precise spectral 
shape, but are similar enough to indicate that the behaviour of the 
source is similar to that in 1997. 



The bright spikes of emission are shown at higher time resolu- 
tion in Fig. 2. The most rapid large rise occurs in Orbit 2, where the 
count rate jumps by by about 4.5 ct s~ 1 from (2.42 — 2.44) x 10 s, 
i.e. 2000 s. Using the spectral model developed in Section 3, this 
corresponds to about 7 x 10 40 erg s~ 2 in the 0.3-2 keV band. An 
even faster event occurs at the peak of Orbit 4 where the count rate 
jumps up and down by 3 ct s _1 in just a few 100 s. This means 
that active regions of the source at that time must be smaller than 
a few 100 light seconds in size. The rate of change of luminos- 
ity is thus approaching the value of 10 42 erg s~ 2 , especially if the 
energy band is extended down to 0.1 keV, which is the highest pre- 
viously recorded, by Brandt et al (1999). That paper discusses the 
lower limit on the radiative efficiency required by the source in 
converting mass to energy to produce such a luminosity gradient 
(Fabian 1979); the efficiency must be high and approaching 50 per 
cent. This is difficult to envisage taking place without invoking non- 
spherical geometry and/or relativistic effects. 

Light curves in both Soft (0.3-1 keV) and Hard (1-4 keV) 
bands are shown in Fig. 3, together with the Soft/Hard ratio and 
background count rate (note that the log of the background rate is 
shown as the values are mostly small). The pronounced large spikes 
in count rate show no strong spectral variation. A trend for the spec- 
trum to be sometimes softer when the source is faint is apparent 
(Fig. 4). Orbit 3 shows considerable hardness ratio variations that 
do not seem to correlate with the count rate of either the source or 
background. 



4 SPECTRAL FITS 

Following the phenomeonological approach of Fabian et al (2009; 
see also Ponti et al 2010), we fit the spectrum over the 0.3-0.4 keV 
and 1.2-2.2 keV bands with a simple absorbed power-law plus 
blackbody model and then show the residuals to that model over the 
full 0.3-9 keV band (Fig. 4). Clear, broad emission residuals cor- 
responding to the iron K and L bands are apparent. A good fit can 
be made over this range with that simple model plus two relativisti- 
cally blurred lines, at 0.92 and 6.7 keV using the LAOR model. The 
blurring parameters are tied between the two lines, yielding an in- 
clination of ~ 60 deg and an emissivity index of 6.6. The equivalent 
width of the lines are 59.6 eV and 1.74 keV. Untying the parame- 
ters leads to a slightly better fit but unphysical energies (1.02 and 
7.06 keV) and a lower inclination (~ 30 deg). 

In practice we do not expect that the emission peaks are 
due to single lines but to line and absorption complexes in the 
reflection spectrum. We have therefore fitted the data with a 
physical model consisting of blackbody, power-law and two 
reflection components, one of high and the other of low ionization, 
similar to the best fitting model for 1H 0707-495. The motivation 
for the two ionization components is to model better a turbulent 
accretion disc. The model used is phabs*zphabs* (blackbody+ 
kdblur* (atable{ref lionx .mod}+atable{ref lionx .mod}) , 
where the relativistic-blurring convolution model KDBLUR acts on 
the ionized reflection model REFLIONX of Ross & Fabian 2005)). 
The results of this fit are shown in Table 2 and Figs. 6 and 7. 

We obtain the best spectral fit if the higher ionization compo- 
nent is replaced by a gaussian line at 0.86 keV (which is then rel- 
ativistically blurred along with the lowly-ionized component). The 
need fro any absorption edge now disappears. This suggests that 
the spectral model we are using is incomplete. It is possible that the 
blackbody component is part of the disc emission and helps heat the 
disc, so increasing the collisional part of the Fe-L emission which, 
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Figure 1. Full band XMM lightcurve (0.3-10 keV) lightcurve of IRAS13224-3809. Bins are 200 s. 




Table 1. Values of model parameters used in emisivity profile determina- 
tion. The absorption component is fixed at the Galactic value. 
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Figure 2. Lightcurves of regions of peak count rate per orbit in the 0.3- 
10 keV band. Bins are 50 s. 



for a temperature of 0.3 keV (the peak of the blackbody) peaks be- 
tween 0.8 and 0.9 keV. Generating the appropriate grids of models 
for this is beyond the scope of the present work. 



4.1 The spin of the black hole 

The spectral fit shown in Fig. 6 requires a steep emissivity profile 
from ~ 1.35r g which, if identified as the innermost stable circular 
orbit (ISCO), means that the black hole is close to maximal spin. 



Component 


Parameter 


Value 


Absorption 
powerlaw 


Galactic Nu cm -2 
Photon index, T 
Norm 


5.3 x 10 20 
r, 7n +0.007 

z - /u -o.oi 
3.86 x Kr 4 


RELCONV 


Inclination, (deg 

Inner Index, q\ 
Outer Index, qi 
Spin, a 


63.8 ±0.4 
2.1 ±0.3 
>9 

i 4 +0.05 
^0.1 

0.988 ±0.001 


blackbody 


temperature, kT keV 
Norm 


0.1 03 ±0.0007 
3.65 x K)- 5 


EXTENDX 


Iron abundance / solar 

Ionization parameter, £,\ 

Ionization parameter, £2 

Normi 

Norm2 

Redshift, z 


> 16 

20.7 ±0.4 

325+^ 
3.0 x 10~ 8 
5.0 x 10~ 6 
4.06 x ur 2 




Z 2 /dof 


955/906 



Fitting the spectrum with the blurring kernel REFLCONV (Dauser 
et al 2010) we find the spin to be 0.988 ± 0.001 (see also Table 1). 

The value of spin is robust to changes in the model with dif- 
ferences in a being less than one per cent in all models tried. A 
small improvement in x 2 is obtained by including a cold iron line 
(rest energy 6.4 keV); a slightly broad one with a — 0.23 keV being 
preferred. Including a small edge at 1.15 keV produces a more sig- 
nificant drop in x 2 ■ This may indicate the presence of a hot wind, 
as found for 1H 0707-495 by Dauser et al (2012), or just be a defi- 
ciency in the model since it occurs where the reflection components 
are dropping steeply (particularly see Fig. 5). The variation of spin 
with inclination is shown in Fig. 8, where no strong correlation be- 
tween these parameters is evident. 

Although the statistical uncertainty on the measured spin is 
well below one per cent, there are larger systematic uncertainties 
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Figure 3. Hard (1-4 keV) and soft band (0.3-1 keV) light curves (top two panels), softness ratio and log (total background rate) (lower two panels). Bins are 
1000 s. 
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Figure 4. Softness ratio plotted versus Soft count rate. 



which have yet to be determined. The most important is perhaps the 
implicit identification of the innermost radius of the reflector with 
the ISCO. Computations suggest that the uncertainty here is small 
and could be less than 0.5r g (Reynolds & Fabian 2008; Shafee et 
al 2008). The work of Schnittmann et al (2012) emphasises that 
emission from matter on plunge orbits is beamed mostly into the 
black hole. We note that the requirement for a low ionization com- 
ponent emphasises that the disc remains dense and thus thin within 
the final gravitational radius. 



4.2 Inferring the position and size of the power-law source 

The break in the emissivity profile at only ~ 2. lr g indicates that the 
power-law source is close to the black hole, within a few gravita- 




Encrgy (keV) 

Figure 5. Ratio of observed spectrum to a model spectrum. The model con- 
sists of a power-law, blackbody and two Laor broad lines which have been 
are fitted to the data. The normalizations of the Laor lines have been set to 
zero before displaying. 



tional radii, and thus must be small and confined within that radius 
(e.g. Wilkins & Fabian 201 1, 2012). 

Confirmation that the source is very close to the black hole 
comes from the reflection fraction. This is the ratio of the reflection 
components to the power-law component, normalized so that unity 
corresponds to a reflector subtending 2n sr. This is not straightfor- 
ward to calculate for a high T source since the REFLIONX model 
does not tabulate the total flux, but only that above 0.1 keV (the 
flux at lower energies is of course included in the computations). 
We assess the reflection fraction by comparing the ratio of the am- 
plitude of the Compton hump around 30 keV of the low ioniza- 
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Figure 6. Full band pn spectrum fitted with double reflectors and a black- 
body component. 
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Figure 7. Components of the best-fitting model shown in Fig. 6. 



tion reflection component to the power law with that predicted by 
the PEXRAV model. The result is a reflection fraction of about 15, 
which is a strong indication of light bending close to the black hole 
(Martocchia & Matt 1996; Miniutti & Fabian 2004). 

The emissivity profile has been determined in more detail by 
fitting the spectrum above 3 keV by the sum of relativistically- 
blurred emission profiles from contiguous radii of the disc (see 
Wilkins & Fabian 2010 for more details; the energy range is re- 
stricted to the Fe-K band as the Fe-L band consists of many over- 
lapping emission lines). The result (Fig. 9, top panel) indicates 
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- IRAS 1 3224-3809, July 201 1 
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Figure 8. 68, 90 and 99 per cent confidence contours for spin and disc 
inclination. 



Figure 9. Top: Emissivity profile obtained by fitting the data in the broad- 
ened Fe-K band with emission from many annuli. Bottom: Integrated pho- 
ton flux as function of radius. 



a triple power-law emissivity profile, which integrates to show 
where the observed photons originate from on the disc (Fig. 9, 
lower panel). 80 per cent of the photons are reflected within about 
2.5r g with the remaining 20 per cent mostly coming from within 
10-20r g . 

Motivated by this, in order to estimate the size and location of 
the major primary X-ray source, we generated a grid of emissivity 
profiles for a range of cylindrical X-ray source regions of varying 
radial extent and at varying heights above the plane of the accretion 
disc using the high speed GPU-based general relativistic ray trac- 
ing code of Wilkins & Fabian 2012. The grid consists of sources 
extending between 1 and 50r g radially, the bases of which are be- 
tween 1.0 and 3.1r g above the plane of the accretion disc. Initially, 
the thickness of the source is set to be 0.5r g . 

These emissivity profiles were then fitted to the profile of the 
Fe-K emission line using a modified version of the KDBLUR con- 
volution model, leading to the constraints shown in Fig. 10. These 
fits imply the source is either radially extended to ~ lr g at a height 
of h ~ 2— 25r g or extended out to around 2 — 3.5r g at a height of 
~1.7r g . 

The integrated emissivity profile of IRAS 13224-3809 is com- 
pared with the normal and low state profiles of 1H 0707-495 in 
Fig. 9. It appears to be sandwiched between the two. This could 
be either due to the source always having two primary emission 
components, one compact as shown in Fig. 10 the second a more 
extended component, or to the source varying in size with time, the 



© 0000 RAS, MNRAS 000, 000-000 



6 A. C. Fabian et al 




Source Radius 



Figure 10. 68, 90 and 99 per cent confidence contours for the radius and 
height (lower edge) of a slab source of thickness 0.5r g . 
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larger source likely being associated with the brighter phases. This 
will be investigated in later work. 



Figure 11. Fractional RMS variability spectrum using 500 s bins. 



5 RAPID VARIABILITY AND THE SOFT LAG 

The fractional RMS variability spectrum, computed according to 
the prescription of Edelson et al (2002), is shown in Fig. 1 1. It re- 
sembles that of many other sources in which reflection is present, 
resembling a combination of variable power-law and reflection 
components. The amplitude of the variability of the power-law 
component needs to be greater than that of the reflection in order 
that the broad Fe-K line appears inerted in this Figure. 

Using the light curves of the four orbits, we compute the 
Fourier phase lag between the hard and soft energy bands, fol- 
lowing the technique described in (Nowak 1999). The background- 
subtracted light curve segments range in length from 8.34 x 10 s 
to 1.24 x 10 5 s with 10 s bins. The soft band is defined from 0.3 
- 1 keV, where the soft-excess dominates the spectrum. The hard 
band, 1.2-5 keV, is dominated by emission from the power law 
continuum. From the Fourier transforms of the hard and soft band 
light curves, 5 and H respectively, we compute their phase differ- 
ence, <j)(f) = arg [(//*£}], where * denotes complex conjugate. We 
convert this to a frequency-dependent time lag, t(f) = <p(f) /2n f. 
Using this sign convention, a negative lag means that the soft band 
light curve lags behind the hard band. 

The results are show in the lag-frequency spectrum in Fig. 1121 
The hard flux lags behind the soft by hundreds of seconds at fre- 
quencies less than ~ 2x 10~ 4 Hz. At frequencies v ~ [3 — 5] x 
100 s. The light-crossing time of 2r g for a mass of 5 x 10 6 M Q is 
~ 50 s, so a total lag of 100 s or so is reasonable. 



6 DISCUSSION 

IRAS 13224-3809 is remarkably similar in overall X-ray behaviour 
to 1H 0707-495. The variability of IRAS 13224-3809 may be the 
most extreme. We shall explore the behaviour of the source as a 
function of time and flux in more detail in later work. 

The X-ray spectra of both sources require high iron abundance 
(Ap e ~ 10 — 20). In recent work, Wang et al. (2012) have presented 
a strong correlation between metallicity, as measured by the Si IV 
O IV ] / C IV ratio, and outflow strength in quasars, as obtained 
via the blueshift and asymmetry index (BAI) of the C IV emission 
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Figure 12. Lag-frequency spectrum for this 500 ks observation. The lag is 
calculated between the soft energy band (0.3 - 1. keV) and the hard band 
(1.2 - 4. keV). We adopt the convention that negative lag mean the soft 
band lags behind the hard band. The most negative lag (at 3.4 X 10~ 4 Hz) 
is -92.1 ±30.7 s. 



line. Their results indicate highly significant super-solar metallic- 
ity (Z/Zq > 5) for quasars with BAI> 0.7. This results indicates 
that metallicity likely plays an important role in the formation and 
acceleration of quasar outflows as expected, for instance, if quasar 
outflows are predominantly line-driven. 

As mentioned above, both IRAS 13224-3809 and 1H 0707- 
495 are characterised by extremely blueshifted C IV emission lines 
with almost no contribution at rest wavelength. Their UV spec- 
tra indicate that BAI> 0.9 in both sources, as shown in Fig. [T3] 
If the metallicity-BAI correlation of Wang et al. (2012) extends 
or saturates above their largest observed BAI (~ 0.76), one in- 
fers that IRAS 13224-3809 and 1H 0707-495 are characterised by 
Z/Zq > 8. A strong indication for super-solar metallicity in both 
sources is consistent with the strong Fell lines in the optical spectra 
and was also inferred by Leighly (2004) via photoionisation mod- 
elling of the UV spectra. 

A ~ 100 s soft lag is detected, which is a direct prediction of 
the reflection modelling used for the source. With the many other 



© 0000 RAS, MNRAS 000, 000-000 



Long XMM observation of IRAS 13224-3809 7 




1600 1650 
Observed Wavelength (A) 



1700 



Figure 13. The C IV emission line profile from the HST-STIS observation 
performed on June 1999 with the G 140L grating is shown in the observed 
frame. Data have been slightly rebinned for visual clarity. The vertical line 
shows the expected wavelength of the C IV emission line for a redshift z = 
0.0658. 



lags now seen, this justifies the reflection spectrum approach. It is 
consistent with the spectral modelling which indicates that the bulk 
of the primary continuum emission source is only a few gravita- 
tional radii in size and distance from the black hole. The spin of the 
black hole is high and close to maximal. This may be the result of 
secular evolution dominating in Narrow-Line Seyfert 1 galaxies, as 
inferred by Orban de Xivry et al (201 1). 
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